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ABSTRACT  OP  THE  DISSERTATION 


Flutter  Control 

with  Unsteady  Aerodynamic  Models 
by 

Shyang  Chang 

Doctor  Of  Philosophy  in  Engineering 
University  of  California,  Los  Angeles,  1984 
Professor  A.V.  Balakrishnan,  Chairman 

This  dissertation  deals  with  a  generic  problem  for 
aircraft*  control  laws  for  flutter  suppression.  Until 
recently,  the  system  frequency  response  was  approximated 
by  rational  functions  so  that  the  finite-dimensional  L- 
Q-R  theory  could  be  applied.  However,  discrepancy 
between  theory  and  practice,  especially  in  transient 
response,  has  led  to  renewed  interest  in  the  problem. 

It  would  appear  that  the  L-Q-R  theory  would  need 
infinite  dimensional  state  space  models.  In  this 
research,  we  first  develop  a  time-domain  model  for 
unsteady  aerodynamic  loads  and  then  couple  it  with  a 
lumped  model  for  the  structural  dynamics.  We  show  that 
the  solution  to  the  resulting  input-output  system, 
characterized  by  integro-differential  equations,  can  be 


endowed  with  a  state  space  which  is  a  reflexive  Banach 
space,  and  the  state  equations  have  a  unique  semigroup 
solution.  We  go  on  to  examine  the  input-output  stability 
for  such  a  system.  We  show  that  input-output  stability 
need  not  imply  stability  of  the  states.  By  a  suitable 
approximation  of  the  Sears  function  near  the  origin,  we 
show  that  the  infinite  dimensional  (L2)  L-Q-R  theory  can 
be  applied.  We  derive  optimal  feedback  control  laws 
ensuring  "weak"  stability  of  the  states,  as  well  as 
input -output  stability. 


CHAPTER  1 


INTRODUCTION 

This  dissertation  is  in  the  general  area  of  modell¬ 
ing  and  control  of  aircraft  flutter  problem.  Unsteady 
aerodynamics  is  also  studied  because  of  independent 
interest. 

NATURE  OF  THE  FLUTTER  PROBLEM 

Aerodynamic  flutter  refers  to  a  subject  that  has 
developed  from  the  earliest  days  of  manned  flight. 

Flutter  is  an  unstable  motion  due  to  the  interaction 
between  structural  vibrations  and  the  aerodynamic  forces 
which  results  in  the  extraction  of  energy  from  the  air. 

It  occupies  a  prime  role  in  current  design  in  the  whole 
spectrum  of  advanced  aircraft »  missiles  and  spacecraft. 

The  field  is  also  one  of  great  inherent  interest  as  a 
scientific  and  technical  discipline. 

The  techniques  required  involve  the  study  of  un¬ 
steady  aerodynamics  for  arbitrary  motions*  structural 
dynamics  due  to  unsteady  loading*  and  aerodynamic  loading 
caused  by  control  surface  motion.  The  primary  design 
goal  is  structural  stability.  Hence*  this  dissertation 
will  focus  on  the  modelling  of  unsteady  aerodynamics* 
coupled  aerostructural  motion  and  flutter  control  systems. 


PREVIOUS  WORKS 


Traditional  methods  of  flutter  analysis  have  pro¬ 
ceeded  in  three  steps i  first  determining  the  vibration 
modes  of  the  structure  without  aerodynamic  forces  pre¬ 
sent  i  then  calculating  the  aerodynamic  forces  on  the  wing 
due  to  simple  harmonic  oscillations  of  the  normal  modes 
as  functions  of  Mach  number*  altitude*  and  reduced  fre¬ 
quency  i  and  finally  calculating  the  flutter  boundary. 

J.W.  Edwards  [l]  and  Edwards  et  al.  [2]  have  ex¬ 
tended  the  unsteady  aerodynamic  theory  from  simple 
harmonic  oscillations  to  arbitrary  motions  using  Laplace 
transform  techniques.  H.  Ashley  et  al.  [3]  and  E.H. 
Dowell  et  al.  [4]  present  similar  formulations  involving 
inverse  Fourier  transforms  to  obtain  impulse  response 

function  airloads.  These  formulations  are  mathematically 
correct*  but  the  calculations  are  cumbersome  and  involve 
functions  available  only  in  tabulated  form.  Hence  few 
examples  of  the  exact  transient  response  of  airfoils  ex¬ 
cited  by  control  surface  motion  have  been  calculated. 

The  first  system  theoretic  formulation  of  the 
problem  is  due  to  Balakrishnan  [18].  Following  his  work. 
Burns  et  al  [6]  have  also  introduced  an  infinite-dim- 
entional  state  space,  although  their  approach  is  based  on 
retarded  functional  differential  equations. 


In  order  to  provide  a  basis  for  the  analysis  of 
aeroelastic  control  schemes,  we  develop  the  exact  tran¬ 
sient  response  of  a  three  degree-of-freedom  airfoil.  The 
development  makes  extensive  use  of  special  time-domain 
functions  derived  from  a  function  studied  by  Kussner 
(Sears,  [5j)» 

Currently  the  active  control  of  flutter  in  flexible 
flight  vehicles  has  gone  beyond  the  frequency  domain 
analysis  in  an  effort  to  apply  linear  quadratic  regulator 
theory  to  the  problem.  However  most  of  the  works  are 
confined  to  finite  dimensional  theory  via  rational  or 
Pade  approximations. 

The  major  problems  with  those  approaches  are  that 

(1)  the  control  obtained  via  finite  dimensional  approx¬ 
imation  is  not  put  back  into  the  original  system i 

(2)  the  singularity  near  the  origin  will  not  be  seen 
after  rational  approximation. 

It  is  apparent  that  the  problem  could  not  be  solved 
without  the  conscientious  introduction  of  infinite  dimen¬ 
sional  state  space.  So  the  first  step  is  to  develop  a 
state  space  model  for  motion  of  an  airfoil  in  two-dimen¬ 
sional  unsteady  flow,  of  an  inviscid,  incompressible 
fluid.  Then  we  show  that  the  solution  exists,  is  unique 
and  depends  continuously  on  the  initial  data.  Of  part- 


icular  importance  in  the  context  of  the  present  work  is 
the  fact  that  we  examine  input-output  statility  and 
L-Q-R  theory  for  such  system  and  thus  may  provide  appro  - 
prite  design  tools  for  flutter  suppression  problems. 

OUTLINE  OF  DISSERTATION 


The  mathematical  model  of  unsteady  aerodynamics  is 
presented  in  Chapter  2.  The  interaction  of  aerodynamics 
with  structure  motion  is  given  in  Chapter  3.  The  com¬ 
plete  equations  of  motion  are  obtained  as  a  set  of 
coupled*  integro-differential  equation.  In  Chapter  4  a 
state  space  model  is  derived  from  the  integro-different¬ 
ial  equation  and  the  well-posedness  of  the  model  is  pro¬ 
ved.  Chapter  5  is  devoted  to  input-output  stability  and 
L-Q-R  theory.  The  final  chapter  presents  the  conclusion 
of  this  research  and  suggestions  for  future  research. 


CHAPTER  2 

MATHEMATICAL  MODEL  OF  UNSTEADY  AERODYNAMICS 


UNSTEADY  AERODYNAMICS 

The  fundamental  equations  of  unsteady  aerodynamics 
follow  from  the  equation  of  state  and  conservations  of 
mass*  momentum  and  energy*  The  derivation  is  presented 
in  numerous  books.  We  follow  the  notation*  nomenclature 
and  development  of  Bisplinghoff »  et  al.  [3]* 


The  exact  equation  satified  by  the  velocity  poten¬ 
tial  is 


where  the  velocity  vector  is  given  by 


q  =  7$  .  (2.2) 

And  the  relationship  between  $  and  a  is  given  by  the  un¬ 
steady  Bernoulli's  equation 


Y  -  1  2  8t  2 

The  boundary  conditions  associated  with  this  pair  of 
partial  differential  equations  are 


DF  _  _3F 
Dt  *  at 


V  F  =  0 


(2.4) 


and  Kutta  condition  for  trailing  edges  of  wings 


The  general  nonlinear  equations  of  potential  flow 
are  difficult  to  solve.  In  order  to  be  able  to  obtain 


solutions  to  these  equations,  it  is  necessary  to  use 
small  perturbation  theory. 

Linearization  is  obtained  by  assuming  that  the  body 
is  thin*  so  that  the  velocity  vector  varies  only  slightly 
from  the  free-stream  velocity.  A  disturbance  velocity 
potential  $  is  defined  such  that 

A 

♦  S  Ux  +  $ 


Then  the  linearlized  partial  differential  equations  for 
unsteady,  compressible  flow  are  (  we  henceforward  drop 
the  *  on  ,  a  ) 


V 


2 


$ 


a 


=  0 


(2.5) 


(2.6) 


subject  to  linearized  boundary  conditions. 


TWO-DIMENSIONAL.  INCOMPRESSIBLE  UNSTEADY  AERODYNAMICS 

As  a  first  step,  we  confine  ourselves  to  incompress¬ 
ible,  inviscid  flow  only.  We  will  derive  a  time-domain 
model  from  the  basic  aerodynamic  equations  and  appropri¬ 
ate  boundary  conditions  (see  Balakrishnan  [?])• 


Consider  a  typical  section  as  shown  in  Fig.  2.1, 
extending  along  the  X-axis  from  -1  to  +1*  with  motion 
entirely  in  the  X-Z  plane.  Then  the  disturbance  velocity 
potential  4>(t,x,z)  is  given  by  the  Laplace  equations 


ifl  ♦  -lIU  ,  0 

ax2  az2 


.  t  >0  f  |  x  I  >  1 


(2.7) 


with  the  boundary  conditions* 


-f^-l  +  -  =  1*1  <  1#  t  > 

a=0  '  a 


(2.8) 


a^tsXjOj  +  y  d*{ttxt0}  _  0  !<x<1+ut  t  >  0  (2.9) 


a»(t.i;o: 


=  0  t  >  0  . 


(2.10) 


The  boundary  conditions  (2.8),  (2.9)  and  (2.10)  are 
the  flow-tangency  condition.  Zero-pressure-discontinuity 
and  Kutta  conditions,  respectively. 

Following  Schwarz  (  see  Bisplinghoff,  et  al.  [3])* 
it  is  easy  to  check  that  the  $(t,x,z)  below  satisfies 
equation  (2.7). 

4>(t,x,z)  - - i-  \  f^iY„(t,s)  tan"1  *  ■  ds  ♦ 


-yw -S'.''.-' ■>  v 


J^Ut  Yw(t.«)  tan’1  J*;  da],  (2.11) 


where  Ya  (t.x)  is  the  circulation  on  the  foil,  and 
Yw  (t,x)  the  circulation  in  the  wake.  They  are  to  be 
determined  from  the  given  boundary  conditions.  The 
integrals  are  to  be  taken  in  the  Cauchy  sense. 


The  nest  step  is  to  represent  4>  in  terms  of  the 
downwash  function  WR  (t.x)  through  the  given  boundary 
conditions. 


“a“  +  _  *  W  (t.x) 

dZ  *-n+»  a 


i  r  r1  i -(t.o 


Moreover, 


flm  1 

ii 


-H  <t.*.0+)  =  \  Y,  (t.*) 


-1  <  X  <  1 


=  I  Yw  (t.x)  1 


<  x  <  1+Ut 


(2.13) 


“§£  (t,x,0+ )  =  0 


x  <  -1,  x  >  1+Ut  (2.14) 


and  that 


(2.15) 


-  -fi  ct.x.o-)  -  -15  <t,x,o ) 


Hence 


4>(t*x,0+)  =  |  j*  Ya  (t.OdC 


,  -1  <  x  <  1 


'  *  ffr- 


(t,C)d£  +  j  Y„  (t,C)dC 


1  <  x  <  1+Ut 


Define 


r(t)  sj^Ya  (t.C)dC 


It  is  clear  from  (2.11)  that 


<t>(t,x,-z)  =  -  4>  (t,x,z,) 


4$  •  -  4r  <*•*•  i*i>. 


In  particular) 


41  <t»x»°  > 


-  "It  (^»x»®  ) 


(2.16) 


(2.17) 


Hence  the  zero-preseure-discontinuity  boundary  condition 
yieldet 


I  [r  (t)  ♦  j*  -Tt  V^>dc  +  UY< 


.(t.x)  =  0  • 


1  <  x  <  1+Ut 


u  Yw  (t.i+)  +  r*  (t)  =  o 


(2.18) 


and  Kutta  condition  yields 


U  Ya  (t,-l)  +  r  (t)  =  0 


Therefore , 


- -  =  y  (t,*l)  »  Y  (t»l+)  and  is  finite. 

U  a 


Prom  (2.18),  Yw  (t.x)  must  satisfyi 


(t.x)  =  P  (  t  -  §  ) 


(2.19) 


Putting  x  =  1, 

Vw  (t.l*)  =  F  (  t 


and  F(t)  =  - 


‘(H) 


.  t  >  o. 


(2.20) 


Therefore , 


Yw  (t.x) 


(2.21) 


v-WSVvV. 


Substitute  (2.21)  into  (2.12) *  we  get 


Wa  (t,x) 


*a(t,o) 
x  -  a 


do 


1  (1+Ut 

*  )l 


r  (t+ffi) 
x  -  a 


* 


or 


Va(t,0) 
X  -  o 


do 


=  -Wa  (t . x ) 


■f 


1 

r 


l*Ut 


r'<t4=2) 

(x  -<j)B 


do 


Wa  (t,x) 


+ 


r’(t-y) 

j  - dy 

J  0  x-l-Uy 


-1  <  x  <  1  .  (2.22 


Since  va  (t,l“)  is  finite*  Sohngen  [8j  proves  that  this 
integral  equation  has  a  unique  solution. 


ya  <t>x)  1  -  hIM  {  fjP  (-*.  <**t> 


+  2"  lo  **)  X-C  dt| 


fVrf 


i  /jjr  w.<t.u 


.jfirr-srrr 


_ if  >I*T  dl 

2”  J  -iV^  *-5. 


This  is  our  basic  integral  equation. 


H  (o.x  )  =  -y=I=uT  •  x-y  ** 


(2.23) 


(2.24) 


and  H  (o)  =  f  H  (o,x  )  dx 


(2.25) 


Then  the  equation  becomes  (  Balakrishnan  and  Edward,  [6]) 


<*•'>  =  ?V¥  { /l.VW 


Wa(t,y) 

x  -  y 


-  5^  Jo  H^»x^  r'tt-hjdn  j 


(2.26) 


where 


r(t) 


=  -2  wa  (t.y)dy  -  j;  j*H(n  )r*  (t-*i)di 


...(2.27) 


Substituting  (2.26),  (2.21)  into  (2.11)  we  obtain  the 
solution  for  the  disturbance  velocity  potential. 


CIRCULATION  ON  THE  AIRFOIL 

With  reference  to  fig.  2.1,  let  h(t)  denote  the 
plunge  coordinate,  <*(t)  the  angle  of  attack,  and  B(t)  the 
flap  deflection.  Then  the  downwash  function  Wft  (t*x)  is 
given  by 


Wa  (t,x) 


tt  s«  <*•*>  +  u  -§*  *,  <*•*>  , 


(2.28) 


where 


•1  «  X  «  c 


(2. 


2"  Jo 


H(<r)r*{t-<r)d(r 


Substitute  (2.30)  into  (2.31),  we  get 


r(t)  =  2irUa(t)  +  (  2oosS+ 2^1-C2  )UP(t  )  +  2trh*  (t ) 


+  (  2yjl-c2  +  (  1-2C  Jcos^c  -  Oyjl-c2  )p*  (t) 


(  rr-2na  )a*(t)  +  )rMt-<r)da 


•••(2*3 


Let 


x( t )  s 


h(t) 

<*(t) 

P(t) 


and 


Z(t ) 


■  f ""  1 


then 


r(t)  = 


[  B  ,  z(t)  j  -*•  J*  (  1-yjEJ  )r* (t-ff)d<r^ 


• • • (2*3 


where 


1'v' 


V.\ 
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V.' 

Ky.  .. 


i 

aaaaa&flJ 


Hence 


-st  tdt 


(2.36 


and 


S',  •-•is 


.  _  1  JS 
z-I  d<r  =  U  e 


K.  <  §  >  “ 


(  §  >  > 


e“Str  r*(*)d <r  = 


U[B,L(Z>s )-  ^p-]eU 

<*,<  *  >  -  <TfTT 


where 


L  (Zis)  = 


Z(t)  dt 


(2.37 


We  will  define  Ci(t)  as  the  inverse  Laplace  transform  of 


1  Ue  U 

8  K.(u>  -  k:<o>. 


Then 


rt  • 

r*  (t)  =  I  C,  ( t- <r)  [B,  z(<r)  I  d<r 


(2.38 


ft 

r(t)  =  J  C^t-ajfB.Z^JJda  -  C2(t)[B,Z(0)] 

•+■  [b,z(o)J  ^  ( 

where 

f  ^ 

c2(t)  =  cx ( s )  ds 

Finally  we  note  the  series  expanision  due  to  Kiissner 
(  Sears,  [5]  )  for  c1(t)  : 

ut<  5 

,  _1  1  3 

0l(t)  =  lAf  H  (Utf2  -  J(Ut)2  +  ife(Dt)2 


In  this  section*  we  will  write  down  the  lift*  pitch 
ing  moment  and  flap  moment  expressions.  The  detailed 
derivations  are  in  Appendices  A*  B*  C. 

1.  P  =  (-P)  {  UIT(t)+  Una*  (t)  -  U  (  C  Vl-c2  -  C0S“1C  ) 

•  p'(t)+  wh"(t)  -  aTra-(t)  -  (|(  Vl-C2)3 

♦  c  cos'^c  -  V  1-c2 )  p"(t)  - 

jf  ^  c3  (t“°  [®*2(Oj  d*  }  (2.42 

where 

c3(t)  =  J*  c1(t-«)  (Uc  -  V  U2c2+2U«]dt  .  (2.43 

2.  Ma  =  (-P)  {  U2(-l-2a/nar  (tj  +  U2  [  (-l-2a)  COS^C 

+  (c-2a)  V^-c2  |P(t)  -  U  (PKZaJirh'U) 

3  3 

+  Ult  (a2-a-j)  tt’{t)  +  U  [  (-|+3ae-2a-2a2-|  ) 

+  (2ac+c-j-2a-a2)  cos^c]  p*(t)  - 

2 

-  (|+a2+a)  tr  h"(t)  -  (j-|-a3-|  )Wa"(t)  - 


2  3  X 

"  a2°“l2“2M  JQ[$  +  Wc^(t-o-) 

X  2 

+  Ua(l-C2(t-<r))j  [B,Z(<r)]d<r  +-|j^[(|+|  ) 

C2(t-a)  +  ac3(t-<r)+  |  C5(t-<r)]  [B,Z(a)]  da} 

. . . (2.44) 

where 

c^(t)  =  c2(t)  +  Cj(t) 

c5(t)  =  |o  ^(t-o)  .  [  (H-Ua)  Vu2o2+2Ua  -  (l+UaJ^io 

...(2.45) 

».  Mp=  (-P)  |U2[(2+c)>/l^c2  -  (2c+l)cos~*c  ]  a(t) 

-  U2fx(c)p(t)  +  U  [  U+c^l-c2  -  (2c+l)cos*1c] 

•  h'(t)  +  U  [(  -  (l+3c)cos~*c 

-  hx(c)]  a*(t)  +  U[^  ”  (^(cj-cf^c))]  P'(t)- 


V  “/  .*  .-  v  v.v 


-  h^fcjh" (t)  +  (ah^cj-h^c))  a " (t )  +  |  (g2(c)  - 

-  cg^cJJpMt)  -  ^  H(Ti)r*(t-ri)dii  +Zfi 


d?  j0  hj 


Oi )  r* (t-njdri 


(2.46) 


where 


fx(c) 


=  ~  [  ^c+lHcos^c)2  -  (1-c2)  -  2  Vl-C2  C08-1c] 


f2(c)  =  ~  {  [(^Ocos^c  +  (c2-c-2)  Vl-o2)3  ] 

.  (  COS^C  +Vl-C2)  +  j(  V  l-c2)^cos“1c  | 


si<°>=  v^"2[  Hs= 


-  -  cos-1®  +  jiuc2 


+  9cf4)]  -  jj-(cos"1c)2(l-2c+2cz) 


g2(c)=  (cos-1c  +^l-c2)(-  J)  [  (1-2c+2c2)cos“1c  - 
^(2c2+9c+4)  ^1-c2  ]  -  |  [(c+^-Jjcos^c  + 

ifc)  ^fl^2  ]  -  [c2fl-c2-  |  cos”1c  + 
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COUPLED  AERO-STRUCTURAL  DYNAMICS 


In  aeroelasticity  the  aircraft  structures  under  con¬ 
sideration  are  perfectly  elastic.  That  is,  the  model 
describing  the  structure  requires  partial  differential 
equations.  As  a  first  step,  we  will  only  treat  a  three- 
degree-freedom  two-dimensional  typical  sections  which  may 
be  regarded  as  representing  the  first  bending  and  torsion 


modest  of  a  three-dimensional  flexible  wing. 


The  equations  of  motion  (  following  the  usual  con 
ventions,  see  Theodorsen  [9] »  and  Edwards  [1])  can  be 
written* 


Bflx  +  B_x  +  K_x  =  zr  +  Gu 


s  m. 


where  the  subscript  s  stands  for  structure*  and 


=  T  1  x . 


x*  vp2”p(c-a) 


i'|«p<c-a)i 


V2  2 

V«wa 


v  2.2 

P, 


U2(-l-2a)tf 

itU2  [(-l-2a)cos“1c 

+(c-2a)^l-c2] 

U2[(2+c) y 1-c2 
-(l+2c)cos“1c] 

-U^c) 

ttU 

-U[c^l-c2- 

cos_1c] 

nU(lf2a)  iru(a2 

0 

-a-|)  u{(-|*3ac-2 

2  c2  2c3 
-2a  -3-— 

Vl^c*-(§-a 

*2a-2ac-c ) 

cos^c } 

[(2<-c)Vl-c2  U[(|+C4^)  U[^“  “ 
(lf2c)cos“1cJ  (f2<c)'c 

3c)cos“1c  *fl(c))] 

(3.2) 


where 


c6(t)=  c1(t-|)H1(e)cU 


r(t  )=  J*  ^(t-DHgUJde  . 


Moreover#  we  can  write  (3»2)  in  a  more  compact  forint 


M 


a.c 


M 


B.c 


=  Jq  M2(t-v)Z(v)dv+  M3(t-v)Z(v)d' 

U[B,Z(0)]1 


(3 


where 


K2<t)=-  ST 


S 


Uc2(t)B* 


M3(t>=-£ 


[^01(t)  4  Uc4(t)  4  «0(l-C2(t)]B* 

-  w*’8* 

[<§  4^)-^e2(t)fac3(t)4jc3(t)]B 

Finally  the  equations  of  motion  can  thus  be  written  as 
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M_X*B_X+K_X  =  Gu+M  X+B  XfK  Xf I  M,(t-6)Z( 6)d6 

8  S  S  EL  8  8  In  2 


ft 

I  M, 

Jo  4 


*  3t|0  M3(t  -6)Z(6)d6+ 


UB*Z(0) j  (  ~"  )  (3.4) 


Z(t)  =  AZ(t)*Hu(t )  +  [B»Z(0)]  fO 


where 


+  0 
0 


jo  <Ms-Ma)“1M2(t-6)Z(6)d6 


t  0 
0 


0  ft 

tt<Ms-Ma)"1]o  M: 


(t-6)Z(6)d6 


A=  0 


(3.5) 


l<w 


I 


Since 


(Ms.Ma)'\(«)Z(°)  =  [B,Z(0)]v 


and 


t  M2(t-v)Z(v)dvsM2(*)2(t)-M2(-)Z(  0) 


M2(t-v)Z(v)dv 


We  have 


Z  =  (A+ 


B#)Z(t  )♦ 


0 

0 

0 


-i/^r^ct-ozMdv 


♦ 


♦Hu(t) 

* 


where 


M2(t)= 


U(c2(t)-1)B* 

[{Jc^tJfUc^tJ+aUU-c^t))^ 


*  25=6 <t)B* 


and 


^  s 


B* 


<§  -a-$)B* 

[ (4c‘fl)cos“1c  -  (2+3c)  ^1-c2]b* 


Now  we  can  rewrite  (3*6)  as 


Z(t)=QZ(t)  <■  j*  K2(t-h)Z0i)dTH>jfJ*  K^(t-n)Z(n )dTi 


♦  Hu(t) 


(3-7) 


where 


k2(«)  =  0 

Moreover  straightforward  calculation  shows  that 

K2(t)  =  0(t"2)  as  t  »  * 

1 

=  0(t“2)  as  t  0 

The  function  K^(«)  is  bounded  and  absolutely  continuous 
on  the  interval  (0,«°).  Thus  let  =  K^(t)  -  K^(«)  and 


carry  out  the  differentiation  in  (3*?)*  we  get 
Z(t)=AZ(t)  *  j*  Q(t-\)Z(x)d\  ♦  rtu(t) 


(3.8) 


where 


-1 

:  5) 

as  t 

~Z) 

as  t 

We  shall  henceforth  work  with  (3*8). 

Finally  in  this  chapter,  we  calculate  the  matrices  in 
equation  (3.8)  by  assuming  the  following  section  para¬ 
meters  * 


U) 

O ( 

=  100 

rad/sec 

a=  -0.4 

Xp  =  0.0125 

=50 

rad/sec 

NO 

• 

O 

II 

O 

r *  =0.00625 

wh 

=300 

rad/sec 

x„=  0.2 

s-° 

=40 

CM 

• 

O 

II 

v 

We  get 

0  1 

H  = 


0 

0 

2243.62 

-43.43 

18977.03 
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CHAPTER  4 
STATE  SPACE  THEORY 


XNTqopugXjQE 

From  the  previous  chapter*  it  was  noted  that  the 
Laplace  transform  of  function  Q(t)  is  not  rational. 

Hence  our  integro-differential  equation  does  not  admit  to 
a  finite  dimensional  state  space  representation.  A  un¬ 
iversal  state  space  formulation  is  due  to  Balakrishnan 
[10] .  a  general  infinite  dimensional  state  space  re¬ 
presentation  for  our  model  is  proposed.  Then  we  show 
that  the  solution  exists,  is  unique  and  depends  con- in 
tinuously  on  the  initial  data.  Moverover  the  controll¬ 
ability  and  observability  issues  are  also  discussed. 

STATE  SPACE  REPRESENTATION 

The  objective  of  this  section  is  to  develop  a  state 
space  representation  for  the  linear  system  described  by 
(3.8)  using  the  techniques  as  indicated  in  [10] . 

Before  a  state  space  representation  can  be  given  for 
(3*8),  some  operators  shall  be  defined  first. 

Let  H^  =  with  the  usual  normi  H2  =Lp-  space  of  6x1 
functions  f(*)  on f 0,»)  with  norm  defined  by 


where  j  =  P  ■*  2 


Prom  now  on  we  will  use  the  Bymbol  IHI  to  denote  any 
one  of  several  norms  when  it  is  clear  from  the  context 
which  norm  is  intended. 

Let  H  denote  the  product  space  H^xHg.  This  will  be 
our  state  space.  Let  A2  denote  the  operator  with  domain 
D(A2)  in  H2  defined  by 

A2f=g  i 

D(A2)=|f«H2  I  f ( • )  is  absolutely  continuous  and 

r<0«H2} 

And  define  operator  B2»  C2  as  follows s 

B2  i  R6  —  H2  B2u=f  i  f(4)  =  Q(4)u 
C2f  =  f(0)  Domain  of  C2=|h«H2  |h(')  is 

continuous  j 

After  defining  operators  A2»  B2#  C2»  it  is  claimed 
that  the  following  set  of  equations  is  a  state  space  re¬ 
presentation  for  (3*8). 


x1(t)=Ax1(t)  ♦  C2x2(t)  -f  Hu(t) 

x2(t) = A2x2(t)  ♦  BgAx^tt)  ♦  B2C2x2(t)  ¥  B2Hu(t) 


These  can  be  rewritten  in  the  forms 


x(t) =  Ax(t)  ♦  Bu(t) 
z(t )  =  Cx(t ) 


(4.1) 


where 


x(t )  =[xx(t)  x2(t)f 


A  = 


B  = 


b2a  a24B2c2 


H 


BzH 


and  C  =[1.0] 


To  verify  this  claim*  see  [lo].  The  main  result  of  this 
section  is 

Theorem  4.1.  The  operator  A  defined  by  (4.1)  generate  a 
CQ-semigroup  on  H^xHg. 


This  theorem  is  an  immediate  consequence  of  the 
technical  Lemmas  (4.1)  -  (4.4)  given  below.  The  proof  of 
Theorem  4.1  will  be  given  after  these  lemmas.  We  will 
concentrate  on  the  "degenerate"  A  (i.e.,  B2=  0)  first. 
Lemma  4.1.  Given  the  "degenerate”  operator 


with  A*  C2t  and  A2  as  dsfined  above .  then  A  is  a  linear 
operator  with  dense  domain. 

Proof.  Straightforward. 

Lemma  4.2.  If  x(*)<D(A2),  then  Hx(0)  II  =  2(llxllp  *►  Ux*llp  ) 

P  2 

Remark t  This  lemma  remains  valid  for  more  general  p. 
Proof.  First  consider  the  restriction  of  x(«)  on  [0,1]. 

x(t)=x(0)  +|  x*(t)  dt  ft  [0,1] 

JO 

or 

x(0)  =  x(t)  -I  x*(t)  dt 

Jo 


Hence 


II x(o) llsllx(t)ll  ♦  I  llx*(t)H  dt 

Jo 

1 

3ix(t)ll  ♦  II x* lip  tq  t  *  [0.1] 


3lx(t)ll  +  II  x  •  ll  p 


nx(0)lips  2p(lix(t)llp  t  llx^) 


Integrate  in  t. 


Ilx(t)llp  dt  «►  llx*llp 

P 

Hence 

tlx  (0)11 p  =  2P  (l!xllp  +  llx*lt  p  ) 
tlx (0  )l!  2  2  (llxllp  *  ltx*llp  ) 


E 


llx(0)ll p  dt 


HE 


I,*mmfl  4.1. 


Q.E.D. 


is  a  closed  linear  operator  on  H. 
Proof*  Given  that 


£  D(A)  S  R6  x  D*A2* 


9 


rx  i 

r„l-i 

A 

rx*~i 

n 

X 

n 

x2 

—* 

x2 

and 

x2 

LxnJ 

•K  « 

LxnJ 

need  to  show  that 


X1' 

«  D(A)  and  A 

X1" 

= 

« - 

>» 

1 _ 

2 

2 

2 

X 

X 

-  « 

7  J 

t 


3? 


It  is  shown  in  Balakrishnan  [12]  that  x2 «  D(A2)  and  A2x2 

2  <v  i  2  1. 

=  y.  We  only  need  to  show  that  Ax  4  =  y.  Accord¬ 

ing  to  Lemma  4.2, 

IlCgX2  -  CgX2!!  =  llx2(0)  -  X2(0)ll 


=  llx2  -  x2ll_  4  ||x*2  -  x*2ll 
n  p  n 


~  l  2 

by  assumption  and  the  result  above.  Hence  Ax*  4  CgX* 
converges  to  Ax1  4  CgX2  and  y^Ax1  4  CgX2. 


Q.E.D. 


Now  we  are  ready  to  show 
Lemma  4.4. 


A  = 


A  C. 

4 

0  A. 


generates  a  C0-semigroup  on  H  =  R^xLp[0,»)^  with  p  ^  2. 
Remark.  The  linear  operator 


A  0 
0  A, 


is  an  infinitesimal  generator  of  a  contraction  CQ- 
semigroup.  C2  is  an  unbounded,  unclosable  operator  by 
itself.  This  lemma  indicates  that 


will  generate  a  semigroup  under  the  perturbation  of  C2. 
Proof.  From  Lemma  4.1. ,  4.3. ,  A  is  linear  closed  and  has 
a  dense  domain.  Moreover*  it  is  easy  to  see  that 

ur( \iA)n  2  xrrnfir  v  x.  >  iiaii 
iir(\iA2)»=^  v  x  >  o 

Set 

R=  R(X|it)  R(\i£)C2R(\»A2) 

0  r(xia2) 

» 

then 


(XI-A)R  =Jxi  -  r  A  0  I  -  TO 

I  Lo  a2J  l  0 

ft*  :>•  t: 

=rxi-*  °  >.r° 

L  o  \i-a2J  Lo 

C2R(  XiA2  )  *| 

0  J 

=  i 

for  every  x  «  D(A), 

R(XI-A)x=  R  (  XI  -  fA  0  >  fO 

C,l) 

=  rUi 


A  0 
0  A. 


,]x-\ 


0  c, 

d 

0  0 


I  R( 
0  I 


\»A  )C, 


0  R(\»A)C, 

I 

0  0 


Therefore  the  resolvent  of 


A  C, 


0  A, 


exists  for  \>HAII  and  it  is  given  by  the  operator  R. 
Next,  we  contend  that  there  exist  some  constants 


M  »» 0  and  w  such  that 


IIR(MA)nlli - for  \>  iu  , 

U-*)n 


n  =1,2,. 


For  n  =  1 , 


*K1  2  =IIR(X  *X)v1  «•  R(uX)C2R(\»A2)v2||J 


'V1 

-v2>- 


Iiv,|i 


+  IIR(\*A2)v2  II  2 

+IIR(\lA)H  ^  II  C*R(\.  |A0) 


v  >>  ii Ail 


Note  here  that 


>-Xt 


-"vP  (}0  •‘wt  A , 


=  llvoll 


2"p 


(4) 


l± 

where 


p-c2 


V  \  >  It  A II 


IIRI,=  rinni 


For  n  =  2, 


2||v1ll  ‘ 


~  (X  HI  All ) 


—V  +  liv2l 

*  H  * 


sAeW*] 


(  x-  ii X n ) 2 


2||V^II 


II  Vo  || 


TT5  4 

All  1  * 


+  1 


(X-IIAII)^  (XHIAII) 
14 


-  H  A»  I 


(  X-llXll 


2  ( X— II  AH  )HIA§ 


2+- 

x  q 


,S_JE. 


I  *  0  S  ~ — rrj 

(  A.-BAU  / 


V  a»-baii 


Proceeding  inductively  we  find  for  general  n» 


M*- 


•V 


: 


tfft 


,v 


F'S> 

m 


k^: 


K.: 


m 


if1 


Vl]  -  2»viii2 
-V2 J  =  (  X — II AH  ) 


II  v2n 


2n  ( \-n  aii  )2n 


r..\2n 


(x-^) 


,  (n  \n~2(  X-1IAII  HliXll**"1)2 

▼  o 


^(n-l)4  -f- 


(x.-iliu)2" 


Lv  2J 


Hence 


HR11  II  = 


V7 


r..\n 


(  X-llAll ) 


for  x  li AH  »  n  =  l,2. 


It  follows  that  A  Is  the  infinitesimal  generator  of 
a  C^-semigroup  S(t),  satisfying  |S(t)|  sV*.""*. 

Q.E.D. 


Lemma  4.5.  A 
A  = 


A  C, 


0  A2+B2C2 


» 


also  generates  a  Cq ■ 


semigroup  on  H  =  R^  x  Lp(0,«o)^  with  j  2:  p  «*  2 


Proofs  First  consider  II  0  0 


i 

I 


0  B2C2RfciA2) 


=  iib2c2r(\»a2)v2ii  siib2ii — ^rnv2ii 


(\q) 


Hence 


JIB,  II 


(\q) 


k 


1  provided  II  B9ll  •«  (xq)q  or 
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v*. 


Next  set 


R=f  R(xiA)fr° 

1r(xia)V 

k*0  \  Lo 

B2C2J  / 

=  R(X»A)  ^1  -J^ 

]r(X»aA 

B2C2J  / 

nB9iiq 

provided  x  ^ -  , 

VII 

c 

V 

CM 

We  contend  that  R  is  the  resolvent  of 
A  C2 

0  VB2  C2 


(XI-A)R 


A 

0 


) 


XX- A 

0 


+  • 


4'-[:  I;  II- 


0  0 
0 


B2C2 


2 

-C2 

XI-A, 


0  0 
0 


B2C2 


R  - 


0  0 
0 


B2C2 
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XI-A  -C2  I 

0  xi-a2J 
0 

\  Lo  b2c2r(xia2 


R(XiA)  R(x.iA)C2R(K»A2) 
0  R(XjA2) 


,¥-{ 


0  0 


0  B2C2R(\|A2) 


(“[•  V 

■  c  :](■  •[ 


o  b2c2r(\»a2, 


0  0 


0  B2C2R(X|A2: 


0  0 


,.](■  -  [ 


0  0 


o  b2c2r(\ia2)J\  Lo  b2c2r(x»a2J 


and  for  every  x  «  D(A) , 


R(\I-A)x 


(  A  c?i  r°  0 

■RiXI'[o  A  H 


0  a2j  lo  b2c2 


■'('"[*  - *[ 


0  0 
0  B,C, 

2  i 


=  £  R(\»A  )( 

k  0 


0  0 
■0  B2c2 


Jr(\iA)  )k  |\I  -  £ 


A  C; 
0  A, 


•7*. 


*  x  . 

Therefore  the  resolvent  of 

A  C2 
0 

exists  for 

llB2Mq 

A  II A II  and  x** — — —  i  2  «  q  *  4  . 

And  it  is  given  by  the  operator  R. 

Finally  we  have  to  show  that  there  exist  sons  constants 
0  and  <*>  such  that 

||  R( Xi  A) US  — —  for  X*.  w  ,  n=  1,2,3* •  •  • 
(X-w) 


Recall  that  in  Lemma 


generates  a  CQ- 


semigroup  S(t)  on  H  with  and  «S(t)H  SV2” e11*11*  . 

For  n-1  , 


u  " 

2 

/  [ 

R 

1 

s 

R(X |A) [ I  - 

u2 

V  l 

R(XiA)  R(\iA)C2R(MA2) 

0  r(x»a2) 


0 

\-i 

Ull 

B2C  2 

R( X»  A) J 

J. 

U2j 

I  0 

0  (I-B^RtxiAg)*1 


...(*) 


Multiplying  it  out,  we  get 


(*) 


R(\»A)u1*R(xiA)C2R(\iA2)(I-B2C2R(X»A2))"1u2 


R( x I A2 ) ( I-B2C2R( xi a2 ) )"Xu2 


2(IR(XiA)u1*24-»R(X»A)C2R(X»A2)(I-B2C2R(x»A2))"1 
•  u2n2)  f  «r(xia2)<i-b2c2r(xia2))"1u2i2 


lU-IIAII)2  1  U-IIAII)2  \  (qx)1/<1  ,  ,,B211  / 

V 


■I  v 


— ) 

7 


II V 


U 


~ f  Hu,!!2  *  - f - J -  )2  „u  „ 

-HAH)2  1  ( A-IIAII)2  Mqx)1/(1“IIBJI  /  2 


"V  ) 

"(q01/q 


H  u2ll 


Note  that 


n  B*iu  t2*{  n  feg  II1  ,1-Vq) 

(<lK)1/q  X_7q)1/,q 

“  (  II B2II  ^3  )  Recall  that  (j)1//q2  for 


x“7f  r  q 


2«*qS4 


(  *1 B2 11  3/4  )  S(  HB2l,1u>«,x  ) 


/  ,  3HB2II  llBpll 
U1 - -±r) - 2- ) 

V  4VF  4VT/ 


vr '  f 

(  for  HB2„.^T  ) 
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ii  h  ii ; 


JT 

\-llSll 


V  X  IIA II  and  X** 


IIB2lq 


or 


HRIi; 


JT 


nB,iiq 

—  V  II  AIM*  — 
II  B„llq  q 


x-iiAn  • 


For  n=  2, 


u. 


rR2(\iA) 


R  (\»A2)[I-B2C2R(\*A2)] 


0 


TV 


r“i1  IK  2IIV2  4  llu2»2 

LuJ  ■  (\-«Xii  J2”  (v-#Xu)2n 


Uq)wqjaJ!|L)2n 

(Uq)1/q  -  «B2i!  )2n 


Hence 


ST  „  ii  B2nq 

- — — - —  v  x *» BAH  ♦  — r —  n  *1,2, 

w  HB9Bq  n  V  q  i 

(X-IAB |— )  , 


It  follows  that  A  C 


0  A2*®2^2 


is  the  infinitesimal 


generator  of  a  CQ-semigroup.  f 

We  are  now  able  to  prove  our  main  theorem. 


Prom  previous  lemmas,  we  have  shown  that 


A  C2 

0  Ag+BgC 


s  an  infinitesimal  generator  of  a  Cft- 


0  A2*B2C2  _  H  B2Hq 

L  j  (IIAKe-  )t 

semigroup  T(t),  satisfying  RT(t)llSf2  e  H 

Observe  that  0  0  is  a  bounded  linear  operator  on 

0  1 

the  reflexive  Banach  space  H,  hence  [  A  C2  1  is 


B2A  A2*B2^2 


is 

. 


the  inifinitesimal  generator  of  a  CQ -semigroup  S(t)  on  H» 


CHAPTER  5 

ACTIVE  CONTROL  OF  FLUTTER  PROBLEM 

INTRODUCTION 

The  vital  role  of  flutter  control  problem  is  played 
by  the  choice  of  aerodynamic  model.  Most  of  the  studies 
up  to  now  were  conducted  by  approximating  the  aerodynamic 
modelling  with  rational  transfer  functions.  As  a  result, 
the  final  model  obtained  is  described  by  linear,  constant- 
coefficient,  ordinary  differential  equations.  The  major 
problems  with  those  approaches  are  as  follows t 

(i)  The  singularity  near  origin  of  the  aerodynamic  model 
will  not  be  seen  after  rational  approximation.  This  im¬ 
plies  that  the  arbitrary  transient  response  predicted 
will  be  less  acurate. 

(ii)  The  control  law  obtained  via  finite  dimensional 
approximations  is  not  put  back  into  the  original  system. 
Hence  the  control  may  not  be  able  to  stabilize  the 
original  system  even  though  it  stabilized  the  approximat¬ 
ing  system.  This  can  be  seen  from  the  following  example. 

AN  EXAMPLE 

Consider  the  following  one-dimensional  version  of 
equation  (3*8) t 


f  t^ 

z(t)  =  z(t)  ♦  I  Q(t-e)z(e)de+  u(t)  (5.] 

Jo 

1 

^  "2 

where  Q(t)  in  our  aero elastic  system  is  of  0(t  )  for 

small  t.  flecall  that  the  rational  approximation  of 

Kiissner  function  is  q  (t )  »  l-O^OOe”0*1^0*-  O.SOOe”1'  in 

engineering  literature. 

2ftn=o.o65=-°-130t  * 


For  this  one  dimensional  example,  we  are  going  to  ill- 
ustrate  that  if  we  take  Q(t)  * e  as  our  corresponding 
rational  aoproximation(  then  the  feedback  control  law 
derived  from  approximate  model  will  not  be  able  to  stab¬ 
ilize  the  original  system. 

First  we  rewrite  the  equation  (5*1)  as  follows t 

z(t)  =  z(t)  w(t)  +  u(t)  (5. 


where 


(o)de 


Differentiating  (5*3) » 

w(t)ss  z(t)  -  w(t) 

»  z(t)  «■  w(t)  ♦  u(t)  -  w(t) 
*  z(t)  ♦  u(t) 


(5- 
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Hence  (5*1)  can  be  written  ast 


x(t) =  Ax(t )  ♦  Bu(t ) 
z(t )  =  Cx(t ) 


(5-4) 


where 


It  is  easy  to  check  that  the  matrix  A  has  one  stable 
eigenvalue 

\  -  i-JT 

1  2 


and  one  unstable  eigenvalue 


Moreover  the  system  (5*4)  is  controllable  and  observable, 
it  is  well  known  that  we  can  find  a  feedback  control 
u(t)=  -kx(t)  such  that  the  system  will  be  stabilized,  if 
we  take 
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w>V 


=  *  .3  x  ♦  .1  +  \s 


8 


This  matrix  has  two  stable  eigenvalues,  hence  the  feed- 


*4 


u 


back  system  is  stable  for  an  approximating  model. 


Now  let's  consider  the  original  system  (5.1)  with 


5(t)=o1(t), 


z(t) =  z(t )  * 


jo$(t-e)z 


z (e )d©-#>  u(t) 


=  z  { t  )  + 


jo$(t-e)z 


z(e)de  -  l.iz(t)  -  o.2w(t) 


&&&&&& ££ 


•  *'•  ^  »'•  A  » 


=  -O.lz(t)  ♦  0.8  J  c1(t-s)z(s)  ds 


Take  Laplace  transform  of  equation  (5*5) » 


8 

■u 


sZ(s)  = -0.1Z(s)*0.8  [sZ(s)-z(0)]  jj-  - 2S - +  z(0 

8  Ko<S>^Ki<§> 


/  0  X 

Z(s)«z(0)f  - jV 

\(s+0.1)(Ko(§)^K1(§))-0.8Ue‘|/ 


z(0 )H(s) 


In  order  to  examine  the  transient  behavior  of  the 
system,  we  take  the  asymptotic  expansions  for  large 
arguments  s  of  the  Modified  Bessel  functions  K0(§)  and 

%<§>• 


where 


_  _a 

Hence,  K0(jj)+Ki(§)  -  V 2-&U  e  U  for  large  s. 

▼  5 

The  denominator  of  the  transfer  function  H(s)  now  takes 
the  following  forms 


J  —  (s+0.1)  -  0.8U 


The  original  feedback  system  is  clearly  unstable.  This 
example  illustrates  that  the  feedback  control  can  stabi 
lize  the  approximate  system  but  not  the  original  system 


input-output  stability 

In  order  to  stabilize  the  original  system,  we  will 
work  with  the  following  problem i 
Given  that 

x(t)  *  Ax(t)  Bu(t) 

z(t) =  Cx(t) 


with 


r 00  f DO 

J[u]=Jq  Jjx1(’t)||P  dt  +  Jol«(^)|P  dt 

we  want  to  find  a  feedback  control  such  that  the  cost 
functional  will  be  finite.  Note  that  the  x^(t)  is  only 
the  finite  dimensional  portion  of  the  state  space. 

We  need  the  following  technical  lemma  in  the  proof 
of  main  theorem  5*1.* 

Lemma  5.1.  (Paley-Wiener)  Let  A  be  a  matrix  with 
entries  consisting  of  positively  supported  delta  distri¬ 
butions  and  L^-functions ,  then  it  is  regular  if  and  only 
if  its  Laplace  transform  satisfying  the  following  con¬ 
ditions 

inf  |det  L(A)|  >0  ,  where  C+*  ]  siRes  £  o( 

8«C  (  ) 

Proofs  See  Hille-Fhillips  [l?] . 

Theorem  5.1.  Given  that  (A,H)  is  controllable ,  then 
there  exists  a  feedback  control  u(t)  such  that  J[u]<«»  . 

Proofs  Prom  chapter  3r  it  is  clear  that  (A»H)  is  con- 


trollable,  then  the  finite  dimensional  Ricatti  equation 

0  has  a  nonnegative  definite,  self- 
adjoint  P1#  such  that  A-HH#P^  is  stable  if  u(t)* 

•x1(t). 

In  order  to  show  that  J[u]  is  finite,  we  need  the 
resolvent  for  A  C2  first. 

BgA  ®2^2*A2 


Recall  that  the  resolvent  for 


A  C. 


is 


A 

R  = 


R(llA) 

0 


R(X|X)C2R(A.iA2)  jl^CgRt^iAg)!"1 
RUtX)  jl-B2C2R( A  »A2)J 


-1 


HB2l,q 

for  s>DJu+  — r — 


Define  Rn*R(;UA) 


R12s  R(AiX)C2R(^iA2)  jl-B2C2R(^lA2) 


-1 


R22 sR(xiA2)  I-B2C2R(A.iA2) 


-1 


Then  A  = 


B2*  B2C2+A2J 


will  have  the  resolvent 


Ai 


R»R£ 

k«0 


0  A  k 
R 


b2a  0 


(I-B2AR12)‘1B2XR11  (l-BjiR12y 


„  UB-llq 

for  A  >ltB -1HI All  ^ II AH  i — f — 
*  <i 


R11‘R12(I“B2AR12)  ®2AR11  Rj^d-BgA  R12^ 


«1  ^ 
I  A 


**22  2^12^  ®2^**11 


R22^  I“B2AR12  ^ 


Now,  let  us  consider  the  following  feedback* 

"x  (t) 

u(t)  =  -kx ( t )  =  -ffi*Plfol  1  *  -H*P,x,  (t) 

L  A  J  x2(t)  1  1 


then 


A-Bk 


B2A  B2C2+A2  B2H6x1 


[«!*6P1, 


X  c. 


B2X  b2c2»a2  B2»r*P1  0 


m*?l  o 


A-HH*P-i 


B2(A-Hll*P1)  B2C2«fA2 


Define  M«  A-HH*P1#  then 


A-Bk 


M 


B2M  B2C2*A2 


•  and  x(t) =  (A-Bk)x(t) 


Note  that  Bk  is  only  a  bounded  perturbation  of  A,  hence 
A-Bk  also  generates  a  semigroup.  In  order  to  show  that 

jj,xl("fc)llP  *  f"lu(t)l  pdt ■«=«,  we  will  examine  the 


behavior  of  R(X»A) 


x^O) 

x2(°) 


instead. 


xx(\)  *  ( B11^B12 ( *“B2^R12 > *1b2®BU )xx (0 ) 


where 


*  N(A iM) t 

R12“  R^»“)C2R<A,A2^I-B2C2R(X,A2)j’1 

Substituting  R^2  into  the  x^(A,)  equation,  we  get 


xl(A)=  [Rll+RnC2R^,A2^  (i’B2C2R^A,a2^)”1, 

•  {^hzyWi]^) 

*  [R11*R11C2R<*»A2>  ((i-b2mr12)(i-b2c2r(a,a2))-1 

B2MRll]xl(0) 
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”[Rll'fRllC2R^A',A2^  I"B2C2R^X,A2^”B2MIb.2+B2MR12B2C2 


R(A.iA2)  “1B2MR11]x1(0) 


“LRllfRUC2R(X,A2 *  jI“B2C2R(X,A2 > -B^5R12( I-B2C2 
•  r(a»a2)|*1b;^r11]x1(o) 

=[R11*R11C2R(*,A2  >  j i“B2C2R( * » A2 ) i A2 ) 
*V*Rll]xl(°) 

=[h11+r11c2r(^,a2)(i*vb2Shii)02r(^,*2) 
t  ( (B2*B2MR22 )C2R(AiA2) )  ♦  •  •  B2®IR^2iJx^ (0 ) 
=[Rn*  Rii^b2**Rii+ ^  b2^b2i®Rii  ^  b2*®r^^+ •  * 


Since  M  is  stable*  the  inverse  Laplace  transfers  of 
Rllxl(°)  will  be  in  Lp[0*»).  Moreover  the  inverse 
Laplace  transform  of  1-(B2+B2MR^)  consists  of  only 
positively  supported  delta  distributions  and  L^-functions 
hence  it  is  regular  if  and  only  if  the  determinant  of 
I-fB-j+BgMR^)  is  bounded  away  from  0  on  C.  by  the  pre- 


A.  A  A  ~ 

vious  lemma*  Since  M  is  stable »  then  ®2*®2*®^11  *8 
bounded  away  from  identity  for  our  aeroelastic  system  and 
the  determinant  condition  is  clearly  satisfied.  The 
second  term  in  (*)  is  again  in  Lp[0,«).  This  shows  that 
J[u]  is  indeed  finite. 

LINEAR  OPTIMAL  CONTROL  THEORY 

In  recent  years*  a  large  amount  of  literature  has 
been  devoted  to  the  L-Q-R  problem  for  infinite  dimension¬ 
al  Hilbert  space.  See,  for  example  [13,  14, .15].  In 
chapter  4,  we  have  shown  that  the  infinitesimal  generator 
A  generates  a  Cg-semigroup  in  a  reflexive  Banach  space. 
The  reason  for  that  is  because  of  the  aerodynamic  energy 
consideration,  i.e.,  the  behavior  of  Q(t)  near  origin  is 
of  0(t~1//2).  If  we  approximate  the  function  §(t)  near 
origin  by  t“1^2+<  ,  then  the  proof  for  Banach  space  will 
go  through  for  Hilbert  space  R^  x  Lg^CO.oo).  As  a  first 
step,  we  shall  in  this  section  work  with  the  state  space 
R6  x  L26[0.oo). 


First  we  recall  the  concepts  of  stability.  Let  T(t) 
denote  a  strongly  continuous  semigroup  over  a  Banach 
space  H.  If  H  is  finite  dimensional  then  Re(®'(A))<  0  is 
equivalent  to  ||T(t)||  -  0  or  for  every  x  in  Hj|T(t)x||  —  0 


The  situation  is  quite  different  in  the  infinite 
dimensional  case.  We  will  gitre  a  few  different  notions 
of  stability. 

Definition  5*1*  A  semigroup  T(t)  is  "exponentially 
stable"  if  ||T(t)||  -  0  as  t  -  «  . 

Definition  5*2.  A  semigroup  T(t)  is  "strongly  stable" 
if  ||T(t )x ||  —  0  as  t  «  for  every  x  in  H. 

Definition  5*3*  A  semigroup  T(t)  is  "weakly  stable"  if 
[T(t)x,y]  — *  0  as  t  —  »  for  every  x,  y  in  H. 

It  is  easy  to  show  that  exponential  stability  im¬ 
plies  strong  stability,  and  strong  stability  implies  weak 
stability.  Note  also  that  T(t)  is  weakly  stable,  so  is 
T*( t ) .  But  this  is  not  true  for  strong  stability,  just 
take  the  left  shift  semigroup  for  example.  However  in 
finite  dimensional  cases,  these  notions  are  equivalent. 
Using  rational  approximation,  J.W.  Edwards (X J  applied  the 

finite  dimensional  L-Q-R  theoty  to  the  flutter  problem. 
However,  this  problem  can  not  be  solved  without  using 
the  infinite  dimensional  state  space  setting  as  the 
reason  is  clear  by  now. 


For  our  aeroelastic  systems,  the  generator  A  will  have 
unstable  eigenvalues  when  the  speed  of  aircraft  is  fast 
enough . 


We  will  consider  first  the  following  L-Q-R  problem. 


x  ssAx(t)  ♦  Bu(t ) 
z ( t )  =  Cx ( t ) 


(5.6) 


where 


A  = 


A  C, 


0  A. 


defined  on  R^  x  t  o,®) 


The  cost  functional  is 


[ujsj^  [x1(t)f  x1(t)J  ♦  u2(t)dt 


where  R  = 


I  0 
0  0 


,  A>0,  which,  for  convenience  hereafter 


will  be  taken  to  be  unity.  Or  we  cam  write  J(ul  as 


j^[x1(t),  x1(t)]  ♦  u2(t)dt 


Before  stating  our  main  theorms  for  this  section, 
we  shall  recall  some  results  of  L-Q-R  theory  in  Hilbert 
spaces  (  see  Balartishnan  [12],  Gibson [15]). 


Let  A»  B  and  R  be  defined  as  in  (5*6)*  with 


Theorem  5.2 

A  the  infinitesimal  generator  of  a  strongly  continuous 
semigroup.  Then  there  exists  a  nonnegative,  self-adjoint 
solution  of  the  steady  state  Ricatti  equation  if  and  only 
if,  for  each  x«H,  there  is  a  control  for  the  initial  time 
s  and  initial  state  x  that  J  [u]<  ®  .  If  P^Cs)  is  such  a 
solution,  then  <Pa)(s)x,x  >H  =  min  Jjs.x^).  The  optimal 
control  u(.)  is  given  by  u(t )= -B*Peox(t)  and  x(t)  =  T(t) 
•x(0),  where  T(t)  is  the  strongly  continuous  semigroup 
generated  by  A  *  A-BB^P^. 

Proof.  We  follow  the  same  notation  as  in  Gibson  Jl4j  . 
First  we  consider  a  sequence  {tn|  where  tn— »  as  n-»  , 
and  try  to  investigate  the  finite  time  L-Q-R  problem  for 
which  tfatn<  For  each  of  these  problems,  we  denote  the 
solution  of  the  Ricatti  integral  equations  by  Pn(.)  and 
the  cost  functional  by  Jn(s,x,.).  An  important  observ¬ 
ation,  which  follows  from  the  fact  that  min  Jn(s,x,u) 
as  <  Pn(s)x,  x>H,  is  that 

Fn<t)*Pm<t)  .  toStstnStm  ....(*) 

Let  x  and  y  be  in  H,  we  have 

2 

<Pn(t)  X,  y>H's  <Pn(s)x,  x>H.<Pn(s)y,  y>H. 


J  (s,x,u  ) *  <  p  (s )x,  x >  —  J (s,x,u)  implies  that 
ri  n  n  w 

sup  j<  P  (s)x,  y  >  H  for  each  pair  x  and  y  in  H. 
n  n 

Then  sup  IIPn(s)l!<«  ®  by  uniform  boundedness  principle, 
n 

Let  {t^}  be  an  increasing  subsequence  of  |tn|.  Then  the 
uniform  boundedness  of  IIPk(s)ll  and  (*)  imply  the  exist¬ 
ence  of  a  unique  self-adjoint  operator  Pw(s)  <  L(H,H) 
such  that  Pk(s)x—  Pw(s)x  strongly  in  H,  x«  H.  Now  we 
have  shown  that  this  is  true  for  a  subsequence.  For  the 
original  sequence  to  hold,  just  use  the  generalized 
Schwarz  inequality.  Q.E.D. 

Theorem  5.3.  Suppose  there  exists  a  control  u(.),  such 

that  for  each  x  c H,  the  cost  functional  J[u]  is  finite, 

and  lim  n x ( t )n  ==  0 *  i.e.,  any  control  drives  the  state  to 
t 

zero  asymptotically.  Then  there  exists  at  most  one  uni¬ 
formly  bounded,  nonnegative  solution  of  the  integral 
Ricatti  equation  on  [0,«). 

Proofs  Let  P( . )  be  such  a  solution  and  define  x(t)  and 
u(t)  as  follows i 

x(t)  =  S(t,s)x  and  u(t)« -B*(t)P(t)x(t) 

Then  limllx(t)ll»  0  and  the  uniform  boundedness  of  ftp(.)n 
t  **« 

imply  lim  <  P(t)x(t) ,  x(t)>H  =  0.  Then 
t  ~*00 


Let  v be  an  admissible  control  for  s  and  x,  then  after 
some  tedious  calculations  (see  Gibson  [15] ),  we  can  show 
that 


J^( s  ,x  *u )  *<P(s)x,  x  > H  =  (  S  ,  X  ,  v  ) 


Thus  u  is  the  optimal  control  for  s  and  x,  and  P(s) 
=  P«o(s),  tn=  s.  , 


Q.E.D. 


Theorem  5.4.  If  f  is  in  L^O,*),  g  is  L2[0,«o)  then 
h  * f  *g  exists  and  h  is  in  L2[0,«>). 

Proof.  See  Dunford  and  Schwarz  [16]. 

Now  we  can  state  a  theorem  for  the  system  (5*6). 

Theorem  5»5»  If  (A,H)  is  controllable  and  (A*,R)  is 
approximately  controllable,  then  the  steady  state  Ricatti 
equation  [Ax,Py]  *[Px,Ay]  ♦  [x.y] -[B*Px,B*Pyl  -  0,  Vx,y«D(A) 
has  a  unique  nonnegative  self-adjoint  solution  # 

and  the  optimal  feedback  control  is  given  by 
u(t)»  -Bn'.xCt) 

Proof.  It  is  a  standard  result  from  the  finite  dimen¬ 
sional  control  theory  that  if  (A,H;  is  controllable,  then 
the  finite  dimensional  Ricatti  equation  A*?^  ♦  PjX  ♦  I 


■/'■V  -V-  V  V  V  F>  V  V  V 


-  ^HH*P1=  0  has  a  nonnegative,  self-adjoint  P^,  such  that 
A  ~  HH*?^  is  stable  if  u*  -H*P1x1.  This  implies  that 


x2(t)»  T(t)x2 


where  T(t)  is  a  left-shift  semigroup. 

Since  (A-?IH*P1)  is  stable  and  x°  is  an  L2-function, 
it  follows  from  theorem  5*^*  that 


11  x,  ( t ) it 2  dt  ♦  f  u2(t)  dt  «*  » 
0  1  JO 


Moreover,  if  (A*,R)  is  approximately  controllable,  then 
x(t)  will  be  weakly  stable  (see  Balakrishnan  [13]).  But 
x2(t)  is  strongly  stable  by  the  fact  that  T(t)  is  a  left- 
shift  semigroup,  even  though  x1(t)  is  only  weakly  stable. 
Yet  in  finite  dimensional  space,  it  is  the  same  as  strong 
stability.  Hence  lix(t)ll  ■++  0  strongly. 


Now  by  theorems  5*2*  and  5*3»»  the  infinite  dimen¬ 
sional  steady  state  Ricatti  equation  has  a  unique  non¬ 
negatives,  self-adjoint  solution  and  that  there  exists  an 

optimal  feedback  control  such  that  u(t)=  -B*P  x(t)  where 

•0 

is  the  nonnegative  unique  solution  of  the  SSRE. 


Q 


D. 
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.£. 


Next,  we  consider  the  same  L-Q-R  problem  for 


x(t)  = Ax(t )  ♦  Bu(t ) 
z(t>=  Cx(  t ) 

where 


A  s 

A 

atm 

«■ 

C2 

.  B  * 

H 

and  C=[l  o] 

V 

®2C2+*2 

b2h 

as  before 
• 

The  same  cost  functional 

J[ujs  j*^fx1(t)#  x1(t)]  dt  ♦  u2(t)dt 

Under  the  same  assumptions  as  before,  we  can  only 
prove  up  to  weak  stabilizability,  but  we  are  only  in¬ 
terested  in  the  finite  dimensional  portions  anyway.  So 
the  finite  dimensional  part  of  the  system  will  be  stabil¬ 
ized  as  before. 
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CHAPTER  6 

SUMMARY,  CONCLUSIONS,  AND  SUGGESTIONS 
FOR  FUTURE  RESEARCH 

Aerodynamic  flutter  refers  to  a  subject  that  has 
developed  from  the  earliest  days  of  manned  flight.  It  is 
an  unstable  motion  due  to  the  interaction  between  struct¬ 
ural  vibrations  and  the  aerodynamic  forces  which  results 
in  the  extraction  of  energy  from  the  air. 

The  vital  role  of  flutter  control  problem  is  played 
by  the  choice  of  aerodynamic  model.  The  design  of  act¬ 
ive  aeroelastic  control  is  greatly  aided  by  the  avail¬ 
ability  of  providing  a  mathematical  model  valid  for 
arbitrary  motions.  In  the  past,  engineers  used  inverse 
Fourier  transform  to  obtain  impulse  response  function 
airloads  for  use  in  convolution  integral  solutions  of  the 
equations  of  motion.  More  common  are  calculations  using 
finite  state,  rational  function  approximations  for  the 
unsteady  aerodynamic  airloads. 

In  this  research,  we  developed  a  time-domain  model 
for  unsteady  aerodynamics.  Hence  the  exact  transient 
response  of  a  three  degree-of-freedom  airfoil  can  be 
obtained. 

Then  we  apply  the  unsteady  aerodynamic  models  to 


foil  as  a  three  degree  -of-freedom  two-dimensional  typical 
sections.  The  upshot  is  that  the  coupled  aero-structural 
dynamics  will  give  rise  to  an  integro-differential 
equation. 

Instead  of  using  rational  function  approximations, 
we  contend  that  the  problem  could  not  be  solved  without 
introducing  infinite  dimensional  state  space.  We  showed 
then  the  solution  exist,  is  unique  and  depends  continue* 
ously  on  the  initial  data  via  semigroup  approach. 

The  principal  advantage  of  using  a  semigroup  form¬ 
ulation  is  that  once  a  system  has  been  shown  to  generate 
a  semigroup,  the  problem  is  well-posed  immediately.  All 
that  is  left  is  the  smoothness  of  the  solution.  This  of 
course  depends  on  the  smooth  properties  placed  on  the 
initial  condition  and  forcing  term. 

Most  of  the  studies  on  flutter  suppressions  are 
based  on  finite  dimensional  L-Q-R  theory.  The  major 
problem  with  those  approaches  is  that  the  control  law  may 
not  be  able  to  stabilize  the  original  system  even  though 
it  stabilized  the  approximating  system.  Hence  we  are 
forced  to  deal  with  the  complex  original  system. 

Pirst,  we  consider  the  input-output  stability 


problem  in  the  reflexive  Banach  space  framework*  then  we 
modify  the  function  Q(t)  near  origin  such  that  $(t)  will 
be  in  1>2  0,»).  Then  we  use  the  machinary  of  L-Q-R 
theory  in  Hilbert  space  to  obtain  the  optimal  feedback 
control  law  by  solving  algebraic  Ricatti  equations  in 
Hilbert  space* 

Obviously  this  research  needs  to  be  continued  in 
various  directions.  One  extension  would  be  replacing  in¬ 
compressible  case  to  compressible  case  of  the  aerodyna¬ 
mics.  Another  extension  can  be  the  replacement  of  lumped 
structural  dynamic  model  by  distributed  parameter  system. 

Finally,  we  can  see  that  all  the  control  theories  in 
Hilbert  space  or  Banach  space  so  far  do  not  quite  fit  our 
needs  for  aeroelastic  systems.  Hence  one  of  the  most 
urgent  extensions  would  be  in  the  control  theoretic 
aspects . 
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APPENDIX  A 


CALCULATION  0?  THE  LIFT 

We  have  in  chapter  2  (see  Balakrishnan,  A.V.  and 
Edwards,  J.W.  f7l): 

P=J  P(x,t)dx 


UYa(t,x)dx* 


y)dy  dx| 


8  f 1 

s(-p){ur(t)-§^l  ^xva(t,x)dx  4.  r* (t)| 


Consider 


X  -  K 


dK.  L.  Pf1  2/T5Tx  H(r,x)dx 
8t 


•  r*  (t-T)dr  * t2  ♦  T2 


Now 


I^T  It 

1-M-  X  -  K 


1 

X-K 


(-h*  (t) 
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U-a)a-(t)  -  u«'(t)dt*  J  V — ZZl 


(-U-c)P  "(t)  -  U (3*  (t )<K 


SS  ♦•nh"  ( t )  -  (aiT-w-wx )  a"  (t)  ♦  Uwcr' (t)  ♦  (C0S~^C 
4Vi-c2)pN(t)  4  (x-c)(cos“1c  +yj-^z I  Log  |.|)9"(t) 


4  up*  (t)  (cos^c  Log  |.|  ) 


where 


Log  |.|  =  log  - 


fcLog  M-- 


\(1-C)(14X)  4 

V(l-C)(l4X)  - 

Vl-C2  1 
x‘c 


l/(lfC)(l-x)_ 
/( l+c)(l-x) 


and  where  we  note  that 


X (l4X ' 


dx*  0*  J  2  x  *og  c  L-c2i 


f1  -  ,  3/2  fl 

J  ^  x(x-c)  log  |.  I  dx  =  £  (1-c2)  ^  x(x-c) 


dx*  |(l4C)  I  j"  ^  X^CI 


Hence 


Next 


Tjs-irh-ft)  ♦  (air)  o"  (t )  -  Uir«»(t)  ♦  (c  COS^C-  yj  1-c2 
+  5  l-c2)3)fj"(t)  ♦  \J(cyjl-cZ  -  cos“1c)  p*  <*t) 


T2*ft  J0  H  (Y.x)dx  r*  (t-v)dY 

-ft  |q  (1  ♦  -  V  U*uc )2-i)r*  (t-c)dc 


Hence 


where 


P»(-P)  UT(t)  •»•  U  ir  a*  ( t )  -  U(c  V^L-c2  -  cos^cjp^t) 
-  aira-(t)  -  (^(Vl-c2)3  ♦  c  cos^c 

-  V  l-c^)p-(t)  -  ft  J0  c3(t-6)[B.Z(6)Jd6 
,(t)-Jo  ^  U262  ♦  2U6^d6 
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APPENDIX  B 


CALCULATION  OF  THE  PITCHING  MOMENT  M. 


We  have* 


V  (-P)(|_i(x-a)UYa(t,x)dx+|E|^(x-a)  J  iYa(t.y)d^ 


*  (-p) 

(|\(x-a)va 

(t.x,dx-i, 

♦  r • 

(t)^  ’ 

> 

) 

f 

J  ^(x-a)Ya(t,x)dx  *-aT(t) 

♦  f  xy_ (t,x)dx 

3-i  a 

fl 

xjdx.f1  kP 

1— X  . 

✓  f1  J  wa(* 

)-i  a< 

tX  ) GX  *  1 

1+x 

nj.rti *- 

•fe£ 

H(e»x)r'  (t-e)d  e^dx 

X-x  X  dx 

l«c  J -I 

1  l*h 

1-h 

W  (t.h)  ft 

Ho 

dl 


-V  U2h2*2U*l)r'  (t-n)dh 


tx  «■  (r(t)  -  r(o);*  J*  c3(t-u)[B,z(u)]du 


T,«  x\\  J-*2*  -h'(t)-(tj-a)«*(t) 

1  J-l  1+x  lJ-l  1-ti 

53sr 


*I^Kc)p,(t)  ■  up(t)l  ^}dx 

nj—^  x  4  U*(t))  -  »(a-l-x)a*(t) 

J-x  14X 

4  U  (COS^C  4  /iflOgH)P(t)  4  |  (C0S_1C4  /l-C2) 


+  (C0S~'LC  4 


log|*|  )(x-c)i  p'  (t)  I  dx 
1-x  > 


r  -▼(*!•  (t)  ♦  Uo'(t ) )  f-naa*  (t)  4  U ( — cos—1c  4  c  ^  1-c2) 
•  p(t)  4  {c  COS”1C  -  |(24C2)  /l-C2^*  (t)  # 

Next,  the  "non-circulatory"  terms  (that  is,  terms  not 
containing  r (t)  )  in« 


»  f1  (x-a)2 

8t  j-l 


Ya(t,x)dx 


*  f1  (x-a)24*(h"(t)  t  Ua'(t))  -  ’r(a-l-x) 

J -1  14X  ' 

’  «"(t)  4  (COS^C  * /US’ log  |.|  )p'(t)  4  (C0S**1C 


J.+X 


/. . 


1-x 


log  I . |  )(x-c) 6" (t) 


2 

(^♦a24a)ir(h"(t)+Ua,(t))  *►  ir(J-|-a3-|-)a"(t)  *  u| 

2  r . . 

♦a*a2)cos_1c  ♦  (g+2a2+^ — 2ac)  ’1-c2|b*  (t)  *|(coe 


+  / l-c2(-*-a24a)  +  f l-c2(4ac2-6a2c-c3-4a-|) 

2 

-  cos“1c(a+|-+^4|4a2c) (t) 


Ka*  (-P)j-aur(t)  ♦  u{r(t)  -  r(o)|  ♦  u  c3(t-n) 

[B,Z(n)Jdn  ♦  u{-nh'(t)  -  Ua(t)w  ♦iraa*(t) 

♦  U  (-cos’1c+c  ^  1-c2)  B(t)  ♦>  (c  cos_1-j(2«»c2)  J  1-c 
2 

8  *  (t ) }  t  r'(t)  -  (^a2*a)tr)hH(t)4Ua*(t)) 

2 

- ’T(§-|-a3-|-)  «”(t)  -  U{(|*a4a2)  cos"1c  ♦  (j*2a2 

2  f  2 

- 2ac)  *  1-c2}b ' (t )  ♦  {(cos^c  (g4|— §-a2c-ac) 


4 /T-c2  ( j+a24§a*^ac2-^a2c-jij-^)  }b"  (t)  ♦ 


♦ 


5t  5 1)^2  H(«.*)to|r,(t-u)du} 


where  the  factor  in  square  brackets  in  the  integrand  in 
the  last  term  is 


jC 


K  f1(x-a)2/^=2-  H(u,x)dx 
J  0  l*x 

*  (kj~~  -  z2  ♦  z  tt)  +  2a(z-  J z2— 1 ) 

^  z-1 

.  a2/_tli  .  1) 

Z-1 

It  is  also  possible  to  split  the  circulatory  and  non- 
circulatory  terms  in  yet  another  way.  Thus  by  not  splitt 
ing 


| 

S' 


A 


•*, 

P 

£ 

£ 


f> 


& 


a)Ya(t,x)dx • -ar(t)  ♦  |  xya(t,x)dx 


but  directly  calculating  the  left  side,  we  have* 

|  ^(x-a)Ya(t,x)dx 

T^*!  i  "  ^x~a)fc*  (t)  ♦  wUa(t)  «-ir  (l-x-a)«*(t) 


84 


4  U(cos“Ac  4/ 


14X 


log  |*|  )p(t)  4  |(cOS“iC  ♦  j  1-C2 

1-x 

4  (cos_1c  4 J — 1--  logM  )(x-c)fp*(t)  idx 
'  1-x  ’ 

*  -(l*2a)irU*(t)  -  (l42a)wh*  (t)  4  2a2**  •  (t)  4  u((-l 
-2a)  cos“1c  4  (e-2a)  A-c2)P(t)  ♦  {(cos'^c  4  J 1-c2 
• (-l-2a)  4  2cos-1c  (ac+|+|*f)  *  ^  J l-c2(l-c2*3ac) } 

•  P’(t) 

Ti, *  ♦  2»r  [  {2-  A2-l  *  a -  alr*(t-e)de. 

JO  z-1 

This  yields  1 

Ma*  (-p){u2(-l-2a)ira(t)  4  U2 1 ( -l-2a ) cos ”1c 
4  (c-2a)  A-c2  }p(t)  -  U(l42a)irh»(t)  4  2Ua2«ra'(t) 

4  U{^(l-c243ac)  4  2cos-1c(ac4^4|4j)  -(l+2a) 

( cos *  A-c2)  }p'  (t)  -  (|♦a24a)JrUo^,  (t) 

-  U{(|4a*a2)cos“1c  4  (j«-2a242£^-2ac) J l-c2fp*(t)- 
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2  2 
-  (cos“1c)(a+g^+^4|^a2c+ac)}p"(t)  ♦ 
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■^-z2«-z  ^z2-l42a(z-  *  z2-1)4  a2 r*(t-0)|d0 


Ka»  (-p)|u2(-l-2a)w®(t)  4  U2 1 ( -l-2a ) cos-1c  4  (c-2a) 
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2  ^  . 
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ujV- 
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♦  a2«l,|r.(t-p)dpl 

/  z-1  1  • 

The  last  two  lines  may  be  written  as 
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APPENDIX  C 

CALCULATION  OF  THE  FLAP  MOMENT  M( 


We  have< 


v\y 


c)P(x,t)dx 
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The  first  term  is  equal  to 


irh*  (t )  -  w(a-l-x)  a*  (t)  ♦  U»*(t) 

Next  let 
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Jc  lex  JO  1-6  1 
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(x-c)\a(t,x)dx 


|  ^(x-cjQliorh'  (t)*n (x-l)a*  (t)*aira*  (t) 

Jc  J  lvx 

HJn^tndx-CF^cJ-cF^cnp'CtJ-UF^cJpU) 

“In  Jo  ^(njrtt-’Ddn 

2(h*  (t)-aa'  (t)*U«(t)  )((1«§)  4.c2-(c4)co3“1c) 

2  r . — ■" 

♦  {  <§-§“)  /l-C2-COs“1c}a*(t)-(F2(c)-cF1(c))p*(t) 


-  UFj  (c)P(t)  -  jj^  ^  ^(Vjr^t-Yjdv 


(x-c)2\a(t,x)dx 


“(irh'  (t)-awer'  (t)*Una(t)  )dx 
lwc 


♦  2  f1(x-c)^CE(x«.l)dxa*(t)  -  (g2(c)-cg1(c)). 


•  8* (t)  -  Ugl(c)6(t)  -  fjp  J  H2U)T*(t-t)dt 

where 

H2(<r)m  §  /“ffl  H(«r,x)dx 

J  c 

W  ^(x-c)2/^-  dx  ji-  d£ 

Jc  1+x  J  c  l-£  *  6 

82(c)'  I  JW/S  to 

JC  1+X  •*  C  I-1!  1 

First  term  is  equal  to 
Zh-^cXh'ft)  -  aa'(t)  ♦  Ua(t)) 

Second  term  is  equal  to 
2h2(c)a'(t) 

Next 

H2(  ff)S  I  ^(X-C)2/^^  H(r,x)dx 
J  c  1+X 


=  (-2^ — s^|(z-c)  (*eos  m2- 


z2-l 


(l4c)  Z-l 


r  /  2 

♦  (z^cHcos^c-  *1-C2)4(  /l-c^-jcos^c-- — 

( -4 ) ( z-c Itan-1— {( z-2c-( z-c )2-i) 

*  z-l  l  2 


'(l+c)  «z- 


r^c  4  (l+2c-|-z) 


Hence  putting  it  altogether  we  havet 


*  U( (2+c)  / l-c2-(2c*l)cos“1c) (h* (t)-aa' (t)4 


4  Ua(t)) 


♦  U  (§♦§“)  ^  1-c2  -  c  cos^c  a*  (t) 


-U(F2(c)-cP1<c) )  (tHTF^cWt) 


X  2 

-  is  J(  hx  ( s. )  r  ( t- o  d&  +ii§sL.  j#  ( t  >  -hx  ( c )  (h"  ( t ) 


g  (c )*C£  (c) 

-ao"(t)4Uo'(t))-h2(c)a"(t)4(-^ — 


4561(0  ]*  H2(u)r*  (t-w)dw 


(-P)  U* j<2*c)'  l-cMac+lJcos'^sJ^tJ-lTf^cMt) 
i-u|(2*c)^  l-c2-(2c+l)cos~1c|h*  (t)4u|(j*>C4^  )  f  1-c2 
(l43c)cos"1c-h1(c)}a*(t)4u{^^-(f2(c)-cf(c)|B»(t) 

>h1(c)h"(t)4(ah1(c)-h2(c)  )a"(t)^(g2(c)-cg1(c))p"(t) 
I?  Jo  »i<wr<t-u>du*§^fej*  H2(u)r(t-u)du|  # 

,(c)  s (  j—-~x  (x-c)2dx 
1  JO  14X 

,(il2|i2£2jcos-lc  _  |(2C2»9c44)  ^ 1-c2 

2(c)«  f  /  *"*  (x»l)(x-c)2dx 
*  JO  14X 

2  j  «— ■ . 

*  (§  ♦j)  COS^C  -  (jg*|^c)  *(l-c2)  . 


APPENDIX  D 
EQUATIONS  OF  MOTION 

The  equations  of  motion  of  typical  section  shown  in 
fig.  2.1  are  derived  from  Lagrange's  equations  (see 
Edwards,  J.W.  [»J  ) 


ft  all  (T-V)  *  3,1 


where  the  Kinetic  energy  is 


z*p(x)  dx 


The  airfoil  deflection  is 


zft=  -h-(x-a)d  -  (x-c)pU(x-c) 


The  potential  energy  V  is 


V.  !<khh2.xao,2.yz)  , 


Thus 


♦  sa  fti*  Sp 


hB+|l| ♦  (c-a)}«  B 


Substituting  T  and  V  into  the  Lagrange's  equation,  we  get 
the  equations  of  motions 


mh  ♦  Saa  ♦  Sp  £  «■  l^h  *  P 


s«*>  *  Ics*  {Ip+sp(c-a,fe  *  *«“•“» 
sph  ♦  |ysp(c-a)ji'+  kp^Mp 

Let 


h' 

P 

X  r 

<* 

.  L= 

Ma 

_P_ 

“0 

m  m 

Then  the  equations  above  can  be  written  as» 

M  x  ♦  B«_x  +■  K_x  *  ~ 
s  s  s  m  . 


(c-a) 


♦Xp(c-a) 


where 


